arXiv:1506.06798vl [cond-mat.str-el] 22 Jun 2015 


Deconfined Quantum Criticality, Scaling Violations, and Classical Loop Models 


Adam Nahum, 1 J. T. Chalker, 2 P. Serna, 2,3 M. Ortuno, 3 and A. M. Somoza 3 

1 Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA 
2 Theoretical Physics, Oxford University, 1 Keble Road, Oxford 0X1 3NP, United Kingdom 
3 Departamento de Fisica - CIOyN, Universidad de Murcia, Murcia 30.071, Spain 

(Dated: June 24, 2015) 

Numerical studies of the transition between Neel and valence bond solid phases in 2D quantum 
antiferromagnets give strong evidence for the remarkable scenario of deconfined criticality, but dis¬ 
play strong violations of finite-size scaling that are not yet understood. We show how to realise the 
universal physics of the Neel-VBS transition in a 3D classical loop model (this model includes the 
subtle interference effect that suppresses hedgehog defects in the Neel order parameter). We use 
the loop model for simulations of unprecedentedly large systems (up to linear size L = 512). Our 
results are compatible with a continuous transition at which both Neel and VBS order parameters 
are critical, and we do not see conventional signs of first order behaviour. However, we show that 
the scaling violations are stronger than previously realised and are incompatible with conventional 
finite-size scaling over the range of sizes studied, even if allowance is made for a weakly or marginally 
irrelevant scaling variable. In particular, different approaches to determining the anomalous dimen¬ 
sions r/y bs and ??Neei yield very different results. The assumption of conventional finite-size scaling 
leads to estimates which drift to negative values at large sizes, in violation of the unitarity bounds. 
In contrast, the decay with distance of critical correlators on scales much smaller than system size is 
consistent with large positive anomalous dimensions. Barring an unexpected reversal in behaviour at 
still larger sizes, this implies that the transition, if continuous, must show unconventional finite-size 
scaling, for example from an additional dangerously irrelevant scaling variable. Another possibility 
is an anomalously weak first order transition. By analysing the renormalisation group flows for the 
non-compact CP n_1 field theory (the n-component Abelian Higgs model) between two and four 
dimensions, we give the simplest scenario by which an anomalously weak first order transition can 
arise without fine-tuning of the Hamiltonian. 

I. INTRODUCTION 

The paradigmatic ‘deconfined’ quantum phase transition 
is that separating the Neel antiferromagnet from the 
columnar valence bond solid (VBS) for a square lattice of 
spin -l/2s. The theoretical arguments of Refs. [IH3j indi¬ 
cate that the Neel-VBS phase transition is described by 
the noncompact CP 1 (NCCP 1 ) model [4], a field theory 
with bosonic spinons z = ( 24 , 22 ) coupled to a noncom¬ 
pact U(l) gauge field: 

£ = |(V - iA)z\ 2 + k(V x A) 2 + pl\z \ 2 + A|z| 4 . (1) 

This theory is defined in three-dimensional Euclidean 
spacetime; the Neel order parameter is proportional to 
z^dz, where d are the Pauli matrices. 

Numerical results for the J-Q model (the Heisenberg 
antiferromagnet supplemented with a four-spin interac¬ 
tion 0) support the validity of this continuum descrip¬ 
tion mm, as does work on the SU(n) generalisation 
of the problem at large n mm- Unfortunately though, 
the existence of a continuous phase transition in both the 
NCCP 1 model and the SU(2) lattice magnets remains 
a vexed question. While simulations of the J-Q model 
are compatible with a direct continuous transition, they 
show strong violations of finite size scaling @ maim mi- 
These persist up to the largest system sizes studied so far, 
and hamper the extraction of meaningful critical expo¬ 
nents pjO]. Additionally, direct numerical studies of the 
lattice NCCP 1 field theory have disagreed as to whether 


the transition is continuous mm, or whether scaling vio¬ 
lations similar to those seen in the lattice magnets should 
be interpreted as the initial stages of runaway flow to a 
first order transition |19j . 

Are the scaling violations seen at the Neel-VBS tran¬ 
sition indeed signs of a first order transition, with an 
anomalously large correlation length naira iisi , or are 
they due to the critical theory possessing a weakly ir¬ 
relevant scaling variable 12112012U, or do they indicate 
something more exotic? This issue remains controver¬ 
sial. Its relevance extends beyond quantum magnets, 
since the critical behaviour of the NCCP 1 model is im¬ 
portant for various other fundamental problems in sta¬ 
tistical mechanics. For example, this field theory is be¬ 
lieved to describe the three-dimensional classical 0(3) 
model when hedgehog defects are disallowed [3], as well 
as the columnar ordering transition in the classical dimer 
model on the cubic lattice [221426] . (In the latter example 
SU(2) symmetry is absent microscopically, but argued to 
emerge at the critical point.) There is also numerical 
evidence that similar scaling violations afflict the SU(3) 
and SU(4) generalisations of the deconfined transition 

Ennui. 

In this paper we introduce a new model which is ideally 
suited for studying the universal features of the Neel-VBS 
transition, and perform simulations on very large systems 
(of linear size up to 512 lattice spacings, and 640 for a few 
selected observables). We verify that the model shows 
the basic features expected from the NCCP 1 field the¬ 
ory (Eq. 0: an apparently continuous direct transition, 
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with emergent U(l) symmetry for rotations of the VBS 
order parameter at the critical point. However, we show 
that scaling violations are even stronger than previously 
appreciated. Conventional finite-size scaling assumptions 
are not obeyed: the data cannot be made to show scaling 
collapse, and quantities that would normally be expected 
to be universal instead drift with system size. The larger 
sizes considered here show that these drifts are stronger 
than the logarithmic form conjectured previously El El- 
In common with Ref. uni, we see a drift in finite-size 
estimates of critical exponents. We show that this is 
more drastic than previously apparent. Estimates of the 
anomalous dimensions of both the Neel and VBS order 
parameters, as extracted from the correlation functions 
G(r) at distances r comparable with the system size (e.g. 
r = L/2) yield negative values at large sizes. Negative 
anomalous dimensions are ruled out for a conformally 
invariant critical point by the unitarity bounds [281 | 2 ijj . 
On the other hand, the decay of G(r) with r for r -C L 
appears consistent with the large positive anomalous di¬ 
mensions suggested for a deconfined critical point. It is 
conceivable that the transition could be continuous, but 
that conventional finite-size scaling could fail as a result 
of a dangerously irrelevant variable m- For example, in 
this scenario correlators G(r) with 1 <r<L would pre¬ 
sumably show the true positive anomalous dimensions, 
while correlators with r of order L would behave anoma¬ 
lously (as in e.g. ((> 4 theory above 4D [30]). The hypothet¬ 
ical dangerously irrelevant variable discussed here should 
not be confused with the much discussed Z 4 anisotropy 
for the VBS order parameter (Sec. IV B) which is dan¬ 
gerously irrelevant in a different sense. 


Therefore - unless there is a reversal of the drift in 
exponents at still larger sizes, which seems unlikely — 
there are two possibilities: either the transition is con¬ 
tinuous with unconventional finite-size scaling behaviour 
(for example as a result of a dangerously irrelevant scal¬ 
ing variable) or it is first order. We will discuss both 
possibilities but cannot rule out either. We do not see 
the conventional signs of a first order transition, such 
as double-peaked probability distributions for the energy 
and other quantities. 

On the other hand, an alternative hypothesis put for¬ 
ward previously — that the scaling violations are due 
simply to a weakly or marginally irrelevant scaling vari¬ 
able mmm — is not supported by our data. We also 
rule out any explanation in terms of unconventional dy¬ 
namic scaling, i.e. deviations from dynamical exponent 
z = 1 : our model has 2 = 1 by construction since it 
is isotropic in three dimensions. This isotropy is also a 
convenient feature from the point of view of simulations. 

Turning to theory, we analyse the topology of the RG 
flows in the NCCP "^ 1 model between two and four di¬ 
mensions, in order to assess the possibility of an anoma¬ 
lously weak first order transition. This analysis unifies 
what is known about this field theory in 4 — e dimen¬ 
sions, in 2 + e dimensions, and at large n, and extends 
previous partial results m ■ Treating n as continuously 


varying, we argue that in 3D there is a universal value 
n* below which the deconfined critical point disappears 
by merging with a tricritical point. The argument does 
not fix the value of n*, which could be greater or smaller 
than two, so does not tell us whether the NCCP 1 model 
has a continuous transition. However, it does have the 
following consequences. If n* happens to be greater than 
two, there is a possible mechanism by which a very weak 
first order transition can appear for a range of n values 
(i.e. a large correlation length can be obtained without 
the need for fine tuning of the Hamiltonian). The argu¬ 
ment also shows that n» is greater than one. Therefore 
the inverted XY transition in the model with n = 1 is not 
analytically connected to the critical point in the large-n 
regime of the NCCP ™ -1 model, contrary to assumptions 
made in previous work. 


This picture for the RG flows also clarifies that the 
usual 2+e expansion of the 0(3) sigma model does not de¬ 
scribe the conventional 0(3) transition in 3D, but rather 
the deconfined critical point (if it exists at n = 2). This 
is natural: hedgehogs are crucial in determining the crit¬ 
ical behaviour of the 0(3) model in 3D [4], and the 2 + e 
expansion presumably fails to account for them. The con¬ 
clusion is also in line with the RG result that the 2 + e 
approach to the 3D O (M) model should fail when M is 
less than a universal value M c , conjectured to be above 
3, as a result of neglecting the topology of the sphere 
[321 . Interestingly, our best estimates for the correlation 
length exponent at the deconfined transition (Sec. |IVE ) 
are close to v = 1 / 2 , smaller than most previous esti¬ 
mates but in good agreement with the 2 + e predictions 
(Sec. |yT4| . 

Returning to lattice models, our numerical strategy is, 
instead of focussing on a simple two-dimensional quan¬ 
tum Hamiltonian, to construct a simple 3D classical 
model that is well adapted to large-scale Monte Carlo 
simulations. While the correspondence with classical lat¬ 
tice models in one dimension higher is a standard tool 
for studying quantum phase transitions, one might at 
first glance think that this tool is not available for de¬ 
confined criticality. This is because deconfinement re¬ 
lies crucially on the fact that the Euclidean action for 
the spins in 2+1 dimensions — unlike the energy func¬ 
tional for a classical spin model in 3 dimensions — con¬ 
tains imaginary terms (Berry phases). The effect of these 
terms is to endow hedgehogs in the Neel order parame¬ 
ter with position-dependent complex fugacities [33ti35j . 
After coarse-graining, this leads to a phase cancellation 
effect that effectively suppresses hedgehogs EH3I- 


Contrary to the naive expectation above, we show 
that the remarkable physics of deconfinement, including 
the suppression of hedgehogs by phase cancellation, is 
present in our 3D classical model. This model is for¬ 
mulated in terms of configurations of loops on a lattice 
and is a variant of the models of Refs. [36] |37|. The 
loop configurations have positive Boltzmann weights, so 
define a conventional classical statistical mechanics prob¬ 
lem. However, the partition function can also be mapped 
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to a lattice field theory for CP" -1 spins, and in this rep¬ 
resentation the Boltzmann weights are not necessarily 
real. We show by a direct calculation that they include 
the complex hedgehog fugacities necessary for deconfined 
criticality. 

The loop model introduced here has qualitative fea¬ 
tures in common with loop ensembles arising in world¬ 
line quantum Monte Carlo techniques for sign-problem 
free Hamiltonians such as the J-Q model [35]. How¬ 
ever, direct simulation of a quantum Hamiltonian leads to 
an ensemble of worldlines in continuous imaginary time, 
whereas the loop model is an isotropic three-dimensional 
lattice model. This is a desirable feature for numerical 
simulations as it fixes an otherwise unknown velocity and 
eliminates a potentially significant source of corrections 
to scaling [35] . The geometric form of the model also mo¬ 
tivates new observables — for example, we find it useful 
to consider some percolation like observables such as the 
number of system-spanning strands and the fractal di¬ 
mension of the loops. 


II. LOOP MODEL 

An astonishing variety of critical phenomena can be 
studied using classical loop gases. The present lattice 
model involves two species (colours) of loops, or n colours 
in the SU(n) generalisation. It has a phase in which 
infinite loops proliferate, and one in which all loops are 
short. The short-loop phase spontaneously breaks lattice 
symmetry because the system must choose between four 
symmetry-related ways to pack the short loops. 

The transfer matrix for loop models of this kind gives a 
correspondence with a 2D quantum magnet on the square 
lattice [37j. The colour of a strand is related to the state 
of the spin (at a given point in Euclidean space-time) 
in the quantum problem. The infinite-loop phase corre¬ 
sponds to the Neel phase: the presence of infinite loops 
is equivalent to the presence of long-range spin correla¬ 
tions. The four degenerate short-loop phases map to the 
four equivalent columnar VBS patterns on the square 
lattice. The schematic correspondence between the loop 
model and the continuum field theory (Eq. |T|) is that the 
two species of loops are worldlines of the two species of 
bosonic spinons (z\ ■ z 2 )- (See Sec. [V] for more detail on 
the continuum limit.) 

The loop model is a modification of those studied 
in Refs. [33 [37], with an additional interaction chosen 
to drive the model through a transition without explic¬ 
itly breaking the symmetry between the four short-loop 
(VBS) states. The model lives on a four-coordinated lat¬ 
tice with cubic symmetry proposed by Cardy [10] ■ This 
‘3D L lattice’ is shown in Fig. [l] (left). Formally it can be 
defined by starting with two interpenetrating cubic lat¬ 
tices, C i and C 2l with lattice spacing 2, displaced from 
each other by (1,1,1): 



FIG. 1. Structure of L lattice. Left: Nodes and links (with 
associated orientations) in the unit cell of the L-lattice. The 
two sublattices of nodes are indicated. Right: One of the four 
equivalent packings of minimal-length loops on this lattice. 


The faces of C i intersect the faces of C 2 along lines: these 
define the links of the L lattice. The L lattice is bipar¬ 
tite; its two sublattices are marked in yellow and black 
in Fig. [I] 

We orient the links of the L lattice such that each node 
has two incoming and two outgoing links, with the two 
incoming links parallel and the two outgoing links paral¬ 
lel, as in Figs. Hi This assignment is unique up to a 
reversal of all orientations. 

Breaking up each L lattice node by pairing the links in 
one of the two ways shown in Fig. [2] gives a completely- 
packed loop configuration. In the simplest case, the par¬ 
tition function is just the equal weight sum over all such 
configurations, with one of n colours assigned to each 
loop (the case of main interest here is n = 2): 

z= E L ( 3 ) 

coloured 
loop configs 


We can of course perform the sum over colours explicitly 
to give Z = ^2 loop n 110 ' 1oo p s . Note that the loops are 

configs 

automatically consistently oriented, since the pairing of 
links at a node is always between an incoming and an 
outgoing link. 

With the above Boltzmann weight, the system is (for 
n < 4) in a phase where infinite loops proliferate [361 ; 371 . 
We wish to add an interaction that drives the system 
into a phase with only short loops. First, consider the 
extreme limit of such a phase, which is a configuration in 
which every loop has the minimal possible length (which 
is six links). There are only four such configurations and 



J L. 

r *r 


(2) FIG. 2. The two possible configurations of a node (ay = ±1). 


C 2 = (2Z + l) 3 . 
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they are related by lattice symmetry. One is shown in 
Rg. a right. 

A general loop configuration is determined by the bi¬ 
nary choice of link pairing at each of the nodes. We 
denote this binary degree of freedom at node r by an 
Ising-like variable ay = ±1. To fix the sign convention 
for o’, let us pick one of the four minimal-length configu¬ 
rations as a reference and declare that all ay are equal to 
+1 in it. The four minimal-length configurations are then 
those in which all the as on the same sublattice (A or B) 
have the same sign. We can define an order parameter of 
the schematic form ip = (a a, <Jb) which distinguishes the 
four short-loop phases. This is the analogue of the VBS 
order parameter in the quantum problem. 

We introduce an interaction between nearest- 
neighbour as on the same sublattice (i.e. between nodes 
of like colour in Fig. [I]): 

z = ex p(-- B )’ ( 4 ) 

coloured 
loop configs 

E = —j\ ay.ay-' + ayoy/ ). (5) 

V (r,r')€.A (r,r')EB / 

(The sum over uncoloured loop configurations is equiva¬ 
lent to a sum over the as.) With this choice, there is a 
direct transition at J c between a phase that has extended 
loops and (<p) = 0, and one that has only short loops and 

( 0 ) ^ o. 

As we would expect from the quantum correspondence 
m , the continuum description of the above model is the 
NCCP" -1 model. In Sec. [v] we show this directly by 
mapping the loop model to a lattice CP ra_1 model and 
coarse graining, paying special attention to the fate of 
hedgehogs. 


III. OVERVIEW OF RESULTS 


We first summarise the salient results of our simulations. 

At the most basic level, they confirm that the loop 
model shows the central features of the deconfined Neel- 
VBS transition, and probes the same universal physics 
as the J-Q [5] and related quantum models. 

We find a direct and apparently continuous transition 
between Neel and VBS phases. Fig. [3] shows the order 
parameters for these phases, for various system sizes L, 
very close to the critical point (details in Sec. IV A). The 
data suggests a single transition. This is confirmed by 
examining finite-size pseudocritical couplings J C (L) de¬ 
termined from various observables (inset to Fig. 0); all 
extrapolate to the same value as L —> oo within error 
bars, so we are confident there is a single transition at 


J c = 0.088501(3). 


( 6 ) 


At small sizes, the estimates of critical exponents are 
compatible with those found in the J -Q model at simi¬ 
lar sizes and in direct simulations of the NCCP 1 model 



FIG. 3. Neel order parameter N (dashed lines) and VBS 
order parameter \(p\ (full lines) as a function of J for several 
system sizes L. Continuous curves are interpolations using 
the multiple histogram method m- Inset: pseudocritical 
couplings J*(L) obtained from various observables [N (blue 
rhombi), \(p\ (down-pointing magenta triangles), crossings of 
4 N & |</?| (grey stars), C (black squares), M a (up-pointing 
green triangles), B v (red circles) and V (turquoise asterisks)]. 


0 0 0 DU HI]. As expected 000 , we see an emer¬ 
gent U(l) symmetry for rotations of the VBS order pa¬ 
rameter ip close to this critical point. (The emergence 
of this U(l) symmetry is equivalent to the noncompact¬ 
ness of the gauge field in the continuum action Eq. [l] 
00 .) Within the VBS phase, the emergent U(l) sym¬ 
metry survives up to a lengthscale £vbs that is paramet¬ 
rically larger than the correlation length £. As for the 
J-Q model 0, the U(l) symmetry is apparent in the 
probability distribution for (p\ see Fig. [4| 

Despite the above features, finite-size scaling proper¬ 
ties at the transition are anomalous in various ways. For 
example, an appropriately-defined stiffness for the Neel 
vector — which would be a universal constant at a con- 



FIG. 4. Probability distribution of (p = (ipi,p> 2 ) for size 
L = 64 and (from left to right) J = 0.0886, J = 0.091, and 
J = 0.096, corresponding respectively to the critical point, 
the VBS phase in the U(l) regime, and the VBS phase in 
the Z4 regime. The crossover from U(l) to Z4 symmetry is 
analysed further in Sec. |IV B| 
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FIG. 5. Spanning number A f s as a function of J, for various L. 
Inset: L-dependence of critical spanning number A4*, defined 
( 1 ) as A4 at J = J c (red dots); ( 2 ) as A f a at the point where the 
slope |dA/” s /dJ| is maximal (black dots). Continuous curves 
are power-law fits, A4* = AL P +B, with exponents p = 0.61(5) 
and p = 0.62(3) respectively. Note log-log scale. 


FIG. 6 . Main Panel: heat capacity C versus J. Upper 
inset: peak height C max with constant Co = 0.9836 subtracted 
to account for background. Fit is power law with exponent 
a/v = 1.52(2). Lower inset: standard deviation in energy at 
J = 0.08850. Fit is with exponent ( a/v — 3)/2 = —0.97(4). 


ventional critical point - increases slowly with system 
size, and the critical exponents estimates also drift as 
the size is increased. Similar features were seen in previ¬ 
ous numerical work on the J-Q model museum, 
but the larger sizes considered here show the scaling vi¬ 
olations are stronger than previously apparent. For a 
detailed picture of the transition, we analyse a variety 
of observables. Violations of finite-size scaling are vis¬ 
ible in almost all quantities and do not decrease as L 
is increased. For this reason, we are unable to fit the 
size dependence of the data near the critical point as¬ 
suming either scaling corrections coming from an irrele¬ 
vant scaling variable (even if it is very weakly irrelevant) 
or logarithmic corrections similar to those considered in 
Refs. Eli. (See Secs. |TVA] [TVCl[rVEp 

Fig.[5]shows the ‘spanning number’ Af s versus the cou¬ 
pling J for various system sizes. A/" s is the average num¬ 
ber of strands which span the system in a given direction, 
and is a measure of the stiffness of the Neel order parame¬ 
ter. Instead of tending to a universal value as dictated by 
standard finite-size scaling, the crossing points A f* drift 
upwards as a power of L: Fig. [5| inset. (See Sec. IV C 
for details.) 


We calculate the correlation length exponent v and the 
anomalous dimensions ?7Neei and ?7vbs using several ob¬ 
servables. In the text, results will be presented with sta¬ 
tistical errors in the last significant digit shown in brack¬ 
ets in the usual way; for reasons that will be apparent, 
we are not generally able to estimate systematic errors. 


Estimates of v obtained from finite-size scaling anal¬ 
yses of different quantities are in reasonable agreement, 
but drift significantly, from v > 0.6 at small sizes to val¬ 
ues around v ~ 0.46 for the largest sizes. In contrast, 
values of v obtained from the variation of the correla¬ 
tion length with distance from the critical point lie in the 
range 0.45 — 0.5 with less dependence on size (Sec. IVE |. 


Strikingly, the behaviour of the correlation functions at 
the critical point suggests different values of the anoma¬ 
lous dimensions ?7Neei and ?7vbs depending on the range 
of r used to extract them. Values obtained from corre¬ 
lation functions at separation L/2 both drift from values 
above 0.2 to values below zero at large sizes. Negative rjs 
violate the unitarity bound 77 > 0 ESI ESSy In contrast, 
there is evidence that behaviour for r <C L is consistent 
with positive values for the anomalous dimensions. We 
note that the use of correlators at separation L /2 to de¬ 
termine r) assumes finite-size scaling, which is a stronger 
assumption than that the continuity of the transition, as 
we will discuss in Sec. IIV A 11 


In view of the above, the transition can only be con¬ 
tinuous if some subtlety invalidates the usual finite-size¬ 
scaling expectations. (Of course in principle there could 
be a drastic change in behaviour at still larger sizes 
L 512, but the data gives no reason to expect this.) 
Therefore it is natural to ask whether the transition is 
first order, with an anomalously large correlation length. 
But while the probability distributions of various quan¬ 
tities show violations of finite size scaling, we do not 
see the standard signs of an incipient first-order tran- 
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sition — double-peaked probability distributions, etc. 
(Secs. IV A4 IVCl IVD). Fig. [bj shows the heat ca¬ 
pacity C, which quantifies the fluctuations of the energy. 
This has a diverging peak at large sizes, as expected for 
a critical point with a positive heat capacity exponent 
(a = 2 — 3^ > 0). The peak only emerges from the 
background at relatively large L. Surprisingly though, 
the peak fits well to a power law after subtracting a con¬ 
stant to account for the background, C max — Co ~ AL a l v 
(Fig-Ill upper inset). This gives v — 0.44, corresponding 
to a divergence ~ L 1 " 52 . This divergence is much slower 
than the L 3 expected asymptotically at a first order tran¬ 
sition. For a more intuitive picture, the lower panel of 
Fig.[6]shows the standard deviation of the energy, divided 
by the volume, at J = 0.08850. For a first order transi¬ 
tion this should saturate to a constant (proportional to 
the square of the difference in energy density between 
the two phases) while here there is no sign of saturation. 
(More details in Sec. IV D ) 

To shed light on these perplexing observations, we 
analyse the topology of the RG flows in the NCCP n_1 
theory in Sec. Eg The topology we find allows for a sce¬ 
nario with an anomalously weak first-order transition for 
a range of n, as a result of a coupling which ‘walks’ (runs 
slowly) in the proximity of a fixed point located at spa¬ 
tial dimension slightly below three. This is one possible 
reconciliation of the above numerical observations. 

A more radical possibility is that the transition is con¬ 
tinuous but disobeys finite-size scaling because of a dan¬ 
gerously irrelevant variable. This was hypothesised for 
the SU(3) and SU(4) cases in Ref. [27j . In this scenario 
we expect scaling violations in correlation functions when 
the separation r of the points is comparable with L, but 
not when r is fixed and L —> oo. In Secs. |IV A 1| [V 11 1 we 
consider this possibility in the light of the data. 

At present we cannot rule out either scenario; we sum 
up the situation in Sec. |VII| 


IV. NUMERICAL RESULTS 



1 10 100 
r 


FIG. 7. Correlators GNeei(r, L) (top) and GvBs( r , L) (bot¬ 
tom) plotted against r, for various L. Straight lines corre¬ 
spond to the estimates ?7Neei = 0.259(6) and ?7vbs = 0.25(3) 
from Eq. [12] 


one may define N to be the probability that a link lies 
on a strand that spans the system in the 2 direction. If 
the transition is second order we expect finite-size scaling 
forms |H1 for <p and N, 


<|£|) = L-( 1+ ^bs)/2 /v , 

N = I / -( 1 +^el)/2j JV (jl/vgjj j 


( 8 ) 

(9) 


A. Neel and VBS order parameters & correlators 

The deconfined transition separates phases that break 
different symmetries. In the VBS (short loop) phase, 
lattice symmetry is broken: this is quantified by the or¬ 
der parameter (p introduced in Sec. [IIJ whose spatially 
uniform part is 



This is normalised so \p\ 2 = 1 for perfect VBS order 
(there are V s ; tes /2 sites on each sublattice). In the Neel 
(infinite loop) phase, SU(2) spin-rotation symmetry is 
broken. In the loop representation, the magnitude of the 
Neel order parameter is the probability N that a given 
link lies on an infinite loop [37]. For a finite system, 


where SJ = J — J c . However attempting a scaling col¬ 
lapse using these forms gives negative rjs and very poor 
collapse. Raw data for the order parameters was shown 
above in Fig. [3] 

1. Correlation functions 

Next we examine the critical two-point correlation func¬ 
tions for p and the Neel vector. In the loop represen¬ 
tation, the Neel correlator is simply the probability that 
two links lie on the same loop m- We denote these cor¬ 
relators Gvbs(u L) and GNeei(r, L ), where r is the sepa¬ 
ration of the points (taken parallel to a coordinate axis) 
and L is the system size. Raw data is show in Fig. [7] 
Conventionally, at J c one would expect 

G(r, L) = L- (1+ ^c{r/L), 


( 10 ) 
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FIG. 8. Neel and VBS correlation functions at separa¬ 
tion r = L/2, plotted against I on a log-log scale. The 
slopes give effective finite size exponent values 77Neei(T) and 
??vbs (L), plotted in left inset. The right inset shows the ratio 
G(L/2)/G(L/4) versus L. 


with different 77 s and different scaling functions c for each 
of the two observables. This would imply a collapse 
when plotting L 1+v G(r, L) against r/L. Here this col¬ 
lapse fails, because the effective values of 77 at small and 
at large distances differ, as we now quantify. 

The full correlation function is relatively complicated 
because it depends on two lengthscales, r and L. There¬ 
fore a standard approach is to examine GNeei(A/2, L) and 
Gvbs(T/2, L) as a function of L. According to (10) these 
scale as and L ~^ 1+VvBS ) respectively. In Fig. [ 8 ] 

(main panel) we plot these correlators against L on a log- 
log scale. The gradual change of slope as a function of L 
indicates a drift in the effective values of 77Neei and tjybs ■ 
The effective values 77^1 (A), rjVBs(L) determined from 
the slope are shown in Fig. [8] (lower inset). Note that for 
large L the estimates for both exponents reach negative 
values. As mentioned above, negative values of 77 VBS or 
7 ?Neei are ruled out for a continuous phase transition gov¬ 
erned by a conformally invariant fixed point, though see 
below. 

Another way to quantify the violation of finite-size 
scaling is via the ratios 


Gvbs(L/2,L) 
GVbs(£/ 4 , L)' 


G N eel(A/2,L) 

GNeel(A/4, L) ’ 


( 11 ) 


which should be universal according to ( 10 ) but instead 


drift significantly with L ; see Fig. [ 8 j upper inset. 

At certain critical points — for example in </> 4 theory 
above 4D — a dangerously irrelevant variable invalidates 
standard finite-size scaling for the correlators. In this 
scenario it may happen that the correlator is conventional 



FIG. 9. Derivatives of Neel and VBS correlators, 
dGNeei(r)/d?' and dGvBsM/dr. Main panels show raw data 
and power-law fits (dashed lines) with ?7Neei = 0.259(6) and 
rjyBS = 0.25(3). Insets show scaling collapse of L 2+V G' versus 
r/L, with the same values for 7/Neel, vbs- 


in the limit L —i 00 (i.e. G(r, 00 ) ~ with 77 > 0), 

but anomalous when r is comparable with L , or even 
with a smaller power of L (421 . Although a priori there 
is no theoretical reason to expect this phenomenon here, 
it suggests examining correlators in the regime r <C L. 
From Fig. [7] it is conceivable that a well-defined power 
law will emerge in the limit L —¥ 00 , although if so the 
convergence in L is rather slow. 

For further insight we examine the derivatives of the 
correlators, G' mel (r) = dGNeeiM/dr and G' VB s( r ) = 
dGvBsW/dr, in Fig. @ Remarkably, these quantities 
show quite clean power law behaviour up to at least 
L ~ 100 with 

77N,ei = 0.259(6) 77VBS =0.25(3), (12) 

and quite good scaling collapse (insets to Fig. [ 9 ]). Rea¬ 
sonable scaling collapse is also obtained for ‘subtracted’ 
correlators, defined as G(r, L) — G(L/ 2, L). Straight lines 
corresponding to the above exponent values are shown in 
Fig- 13 for comparison with the raw data. 

These results indicate that the strongest effect of the 
scaling violations is on the zero-modes of the fields, which 
have anomalously large fluctuations. This is suggestive 
because in </> 4 theory above 4D, the violation of finite- 
size scaling is caused by anomalously large fluctuations of 
the field’s zero-mode [SOI- The contribution of this mode 
to the two-point function depends on L but not on r, 
so scaling can be repaired by differentiation/subtraction. 
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FIG. 10. Susceptibilities for Neel and VBS order parameters. 
Inset shows peak heights as a function of L. 


The fact that this works perfectly in </> 4 theory is expected 
to be a special feature of the fixed point being free (and of 
the choice of correlator). Nevertheless, the good scaling 
of and Gy BS is striking given the strong violation 

of scaling for the correlators themselves (Fig. [8]), and may 
possibly indicate that a dangerously irrelevant variable is 
playing a role in the scaling violations. 



FIG. 11. Uy bs versus J for several system sizes ranging from 
32 to 512. Inset: maximum value of dVvBs/dJ versus L. 


3. Susceptibilities 

To compare with the estimates above, we calculate ?7Neei 
and ?7vbs from the Neel m and VBS 011 susceptibili¬ 
ties. These are shown in Fig. [To] According to finite-size 
scaling the peaks should diverge as L 2_?? N4al ’ VBS . The in¬ 
sets show log-log plots of the peak heights against L. The 
slopes indicate a downwards drift from ?7vbs = 0.164(13) 
to 77 vbs = —0.35(10), and from r/Neei = 0.355(9) to 
??Neel = 0.126(3). 


2. Fractal structure of loops 


The geometrical interpretation of the anomalous dimen¬ 
sion j^Neei is in terms of the fractal dimension of the loops, 
which according to conventional scaling relations is given 
by df = (5 — 77Neei)/2, and determines the power-law 
relation between the root mean square end-to-end dis¬ 
tance R of a strand and its length (see Ref. m for de¬ 
tails). This again gives a large positive r/Neeh in con¬ 
trast to the drift towards negative values seen in the es¬ 
timate from G]sfe e i(A/2). The simplest fit, taking strands 
with R < 100 to minimise effects of finite R/L , gives 
??Neei = 0.42(6) (data not shown). We note that this is 
considerably larger than Eq. [l2j However, attempting to 
include finite R corrections in the fit gives smaller values 
in the range 0.25 < r/Neei 0.42 [43] . 


4. Binder cumulant 

Fig. [TT] shows the Binder cumulant for the VBS order 
parameter, defined as 


Gvbs = 2 — 


(M 4 ) 

(M 2 ) 2 ' 


(13) 


At a first order transition there should be a dip in I/vbs 
which diverges with the system size [45]. In our case 
there is no sign of this. 


In the inset to Fig. 0 we plot the maximum value 
of the slope dt/vBs/dJ for each L. For a second order 
transition this diverges as L 1 /" at the critical point. From 
the inset we see that there is different behaviour for small 
and large system sizes, giving v = 0.62(1) f or size s L < 64 
and v = 0.476(18) for L > 256. (See Sec. IV E for other 
estimates of v.) 
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FIG. 12. Data for (cos(40)} versus J in several system sizes. 
Inset: (cos(40)} as a function of the scaling variable lA 1 ' 4 (J — 

Jc)F 2 (J-J c ). 


B. Emergent symmetries 


The deconfined criticality scenario assumes that the lead¬ 
ing operator which breaks the symmetry for rotations of 
if from U(l) down to Z 4 is dangerously irrelevant: irrel¬ 
evant at the critical point, but relevant within the VBS 
phase [2] ■ This leads to the prediction of a crossover be¬ 
tween U(l) and Z 4 symmetry within the VBS phase, on a 
lengthscale £vbs which is parametrically larger than the 
correlation length: £vbs ~ £ 1 + M/ 3 ; where 2/4 < 0 is the 
RG eigenvalue of the fourfold anisotropy [46] . This has 
been confirmed in the J-Q model [3112- 

Fig. [4] gave visual evidence for the emergent U(l) sym¬ 
metry in the loop model. A quantitative measure of Z 4 
anisotropy is (cos40), where (p = |<^|(cos0,sin0). Fig. 12 
shows data for sizes up to L = 200. Ignoring scaling 
violations, the anisotropy should behave as [45 ) 


(cos 4 6) = f 5J F 2 


(14) 


where F 2 (x) = 1 + ax + bx 2 takes into account nonlin¬ 
ear dependence of the scaling variable on J (needed here 
because of the larger range of 6 J studied for this observ¬ 
able) and V 4 = v(l + |y 4 |/ 3 ). The inset to Fig. 12 shows 
the attempted scaling collapse using Eq. [14] The expo¬ 
nent V 4 = 1.09(6) is obtained from the fit. This confirms 
the irrelevance of fourfold anisotropy to the behaviour at 
the transition and to the explanation of the scaling vio¬ 
lations. The corresponding value of | j/ 4 1 is dependent on 
the assumed value of v, but is considerably larger than 
the estimate in Ref. [Tj for a variant of the J-Q model. 

The closeness of the finite-size effective values of ?7vbs 
and TjNeei in Fig.[ 8 ]and Eq. [X2|makes it tempting to specu¬ 
late about a much larger emergent symmetry — an 0(5) 


symmetry relating the Neel and VBS vectors. This can 
be incorporated into an alternative field theory for the de¬ 
confined critical point [45] which was argued to be equiv¬ 
alent to Eq. [T] .Ml. This symmetry enhancement would 
be analogous to the emergent 0(4) symmetry of the ID 
spin-1/2 chain, which relates the spin-Peierls order pa¬ 
rameter and the Neel vector [50]. In the future it would 
be interesting to check explicitly for 0(5) symmetry. 


C. Neel stiffness and spanning strands 

A useful observable is the spanning number, A f s , defined 
as the number of strands that span the system in (say) 
the z direction. Its mean value, (A/), may be taken as 
a definition of the stiffness of the Neel order parameter. 
[52l At a conventional critical point, (A f s ) has scaling 
dimension zero and the scaling form 

(N s ) =h(L l / v 5J). (15) 

Therefore (N s ) should be a universal constant at a crit¬ 
ical point, modulo corrections due to irrelevant scaling 
variables: plots of (N s ) versus J for different L should 
cross at J c . In the VBS phase (N s ) tends to zero expo¬ 
nentially in L, and in the Neel ordered phase it grows as 
L. 

The mean spanning number was shown in Fig. [5] Con¬ 
trary to the above expectation, (N s ) appears to diverge 
slowly with system size at the critical point. This is man¬ 
ifested in the upwards drift of the crossing points in the 
main panel. In the inset we show pseudocritical values 
N*(L) defined in two different ways. The data cannot 
be fitted with conventional scaling corrections from an 
irrelevant variable, i.e. N* = A/” s crlt — A L v with negative 
y: attempting such a fit leads to a positive (relevant) y. 

Previous work on the SU(2) J-Q model found a drift in 
a closely related winding number, and proposed that this 
indicated logarithmic corrections to scaling |9]. Similar 
drifts were found for the SU(3) and SU(4) J-Q models 
E2> fitting slightly better to a power law than a loga¬ 
rithm. The larger sizes considered here for the SU(2) 
case show that the divergence is certainly stronger than 
logarithmic. On attempting to represent it by a pure 
power law N* ~ AL a we find that the exponent drifts 
upwards, but a power law plus constant 

N* ~ AL a + B, a ~ 0.6 (16) 

fits the results for all L. This divergence is of course still 
slower than the linear behaviour expected asymptotically 
at a first order transition. 


1. Drift of critical probability distribution 

In addition to the mean ( N s ), we examine the full prob¬ 
ability distribution of N s . Let P ^ be the probability that 
N s is equal to 2 k, meaning that k oriented strands span 
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FIG. 14. Parameter-free scaling: Pi(J) as a function of the 
FIG. 13. Probability P k that the spanning number A f s equals average number Af s {J) of spaning curves. Inset: size depen- 
2k, for a range of L. In the main panel, J is tuned so that dence of the maximum of Pi(J) on a double logarithmic scale. 

Pq = 0.3 for each L. In the inset J is fixed, J = 0.08850. 


the system in a specified direction and k in the reverse 
direction. This again has scaling dimension zero, so con¬ 
ventionally we would expect the scaling form 


P k = g k (L 1 / 1 '6J). 


(17) 


By contrast, at a first order transition, where the short 
loop and infinite loop phases coexist, P k would have a 
peak at k = 0 from the short loop phase and a peak at 
k oc L from the infinite loop phase. 

The distribution P k obtained numerically is shown in 


Fig. 13 for various L. To compare different sizes, we 
tune J for each L so that Po = 0.3 (using the Ferrenberg 
method [53]). For comparison, the inset shows the dis¬ 
tribution at fixed J = 0.0885, very close to the critical 
point. Contrary to Eq. [IT] the data shows no sign of 
tending to a universal distribution. On the other hand, 
neither do we see a double-peaked structure developing. 

The scaling form Eq. [17] would imply that P k is a uni¬ 
versal function of (A/)) for each k. [Explicitly, P k = 
(g k ° h~ x ) ((A4)).] Therefore a plot of e.g. Pi against 
(A4) would show scaling collapse without the need to ad¬ 
just any parameters. (See Ref. [57] for successful exam¬ 
ples of such scaling collapse for the compact CP 1 [i.e. 
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0(3)] and compact CP 2 models.) It is clear from Fig. 
that such a collapse will not work here. Fig. [14] shows t 
for Pi (interpolating curves are obtained with the Fer¬ 
renberg method |53] ). The dramatic failure to collapse is 
quantified in the inset, which shows the maximum value 
of Pi as a function of L. At a conventional critical point 
this would reach a finite constant as L —> oo, while at 
a first order transition it should tend to zero. There is 
no evidence for saturation over this size range, though it 
cannot be ruled out. 

We note that the maximum of Pi decreases below the 


universal value for the 0(3) universality class [37] , which 
is close to 0.4 (a stiffness in the J-Q model also drifts 
beyond the 0(3) value [T?]). This is further confirmation 
that we are dealing with a direct transition rather than 
two separate transitions that are too close to be resolved. 


D. Energy distribution 


The behaviour of the heat capacity, proportional to the 
variance in the energy E, has been discussed in connec¬ 
tion with Fig. 6. Further information is contained in 
the full probability distribution for E (which is defined 
in Eq. |4|. The top panels of Fig. 15 show how this dis¬ 


tribution evolves as J is varied. The bottom panel shows 
the distribution at J ~ J c for various L. We do not see 
a double-peaked distribution. The width of the critical 
distribution also decreases with increasing system size, 
contrary to the expectation for a first order transition 
(Fig-|6| lower inset). 

The Binder parameter 


v =(E 4 )/(E 2 ) 2 - 1, (18) 

plotted in the inset to Fig. [16] is an alternative quantity 
for analysis. The data shows a peak near the transition, 
with a height V max that decreases with L (it is neces¬ 
sary to use the Ferrenberg interpolation method [53] for 
accurate estimates of Vm ax ). At a first-order transition 
Vmax should saturate to a constant as a result of the 
double-peaked energy distribution, while at a continuous 
transition with a > 0 the peak height should tend to zero 
as I4iax ~ L 2 '" -6 . Here, a direct estimate of v using the 
slope gives a value that drifts from v ~ 0.621(5) for sys¬ 
tem sizes 32 < L < 64 to v ~ 0.481(12) for L > 256. 
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FIG. 15. First three panels starting at top: energy distribu¬ 
tion functions for system sizes L = 320, 400 and 512 respec¬ 
tively, for various J. Bottom panel: energy distribution for 
J = 0.08850 for various L. 


However, in addition to the peak, V has a large back¬ 
ground contribution scaling as L~ 3 (see inset). It is nat¬ 
ural to subtract such a correction. This allows I4iax to 
be fitted to the the power law form corresponding to 
v = 0.468(6) — see Fig. 16 main panel. 



FIG. 16. Vlnax is shown as a function of system size on a 
double logarithmic scale with the background contribution 
substracted. Inset: raw data for V plotted as a function of J 
on a semi-logarithmic scale for several values of L. 



FIG. 17. Size dependence of the effective correlation length 
exponent u e g obtained from finite-size scaling of various quan¬ 
tities. 


E. Correlation length exponent 

In this section we discuss two different approaches to de¬ 
termining the correlation length exponent zz, one rely¬ 
ing on finite size scaling and one not (see Sec. IV A for 
anomalous dimensions r/Neei, Wbs)- 

First, we obtain estimates using the standard finite- 
size scaling forms for various observables. For exam¬ 
ple the Binder cumulant for the VBS order parame¬ 


ter (p would naively scale as Uy bs = fu{L 1 / u 5J), so 
that the maximum value of dCA/Bs/dJ should grow as 
L. Therefore we can define an effective exponent via 
i/ e ff(L) -1 = dlog([I/y BS ] max )/dlogL. We calculate such 
numerical derivatives using four consecutive system sizes. 

Fig. 17 shows the resulting estimates v e s{L) from 
Uy BS> from the probability Pq of having no spanning 
strands, and the order parameters N and \(p\. v c s(L) 
drifts from large values, around 0.62 (in accordance with 
previous studies urn usd to values around 0.46. The 

































V. FIELD THEORY FOR LOOP MODEL; 
HEDGEHOG FUGACITIES 
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FIG. 18. Estimates of the correlation length obtained from 
(1) the exponential decay of the spanning number M a in the 
VBS phase, J > J c (red points); (2) the coefficient of the 
linear growth of M s in the Neel phase, J < J c (black points) 
Note: the overall normalisation is arbitrary in the latter case. 


latter is in agreement with the estimate from the heat 
capacity with the background subtracted (Sec. III). De¬ 
spite the drift, is e g(L) remains above the unitarity bound 
2/5. A drift in v was also identified in Ref. |10j . 


The above estimates all rely on finite-size scaling forms 
for observables on the scale of the system size. To avoid 
the assumption of conventional finite-size scaling, we also 
estimate v directly from the correlation length £ in the 
regime where £ -C L. For values of J in the VBS 
(short loop) phase, we determine £ by fitting the In¬ 
dependence of the spanning number to the expected form 
M s oc In the Neel phase, the spanning num¬ 

ber is expected to grow as M s ~ AL/£(J). This allows 
us to determine £( J) up to the overall constant A. 


The results are shown in Fig. 18 The power-law fits 
shown give v = 0.477(4) for the data in the Neel phase 
and v = 0.503(9) for the data in the VBS phase. These 
values are close to the estimates in Fig. 17 at the largest 
sizes. But it is remarkable that here this behaviour sets in 
at much smaller lengthscales. For the VBS phase (where 
we can determine £(</) without the complication of an 
overall constant), the above exponent fits the data well 
starting from scales as small as £ ~ 15. 


We believe that if the transition is indeed continuous, 
the correlation length exponent is close to v = 0.5, so 
considerably smaller than most earlier estimates from J— 
Q and NCCP 1 models. This value is close to 2 + e ex¬ 
pansion results for the CP 1 nonlinear sigma model which 
should apply to the deconfined criti cal p oint assuming 
the transition is continuous (see Sec. VI 4). 


Models for completely-packed, oriented loops can be 
mapped to lattice CP” -1 models with an unconventional 
but simple form El S3- At first sight the continuum 
limit of these models is simply the (compact) CP n-1 
sigma model. Here we discuss this continuum limit 
in more detail, and show that hedgehog defects in the 
CP 71 ' 1 spin configuration contribute imaginary terms to 
the action that are analogous to the Berry phases in the 
Euclidean action for the 2D quantum Heisenberg model 
[351 [5Tj . These terms are crucial for the present model. 
By the reasoning of Refs. w, they change the effective 
continuum description from the compact CP" -1 model 
to NCCP" -1 . (By contrast, the imaginary terms were 
unimportant for the transitions discussed in Refs. 0HE3 
as a result of the lower lattice symmetry there.) 

The quantities of interest are determined by symmetry, 
so it is enough to consider the case J = 0, where the 
lattice held theory for the loop model is simplest. The 
CP" -1 spins are placed on the links l of the lattice. They 
are complex vectors z t = (zj ,..., z"), with fixed length 
|z| 2 = 1 and the gauge redundancy zi ~ e llpi zi. In a 
loose notation where the incoming links at a given node 
are denoted i and i! and the outgoing links o and o', the 
partition function is 


Z = Tr J] 


nodes 




(19) 

Here ‘Tr’ is the integral over the zs. Under a gauge trans¬ 
formation of zi, the terms for the two nodes adjacent to 
l pick up opposite phases, so the Boltzmann weight is in¬ 
variant. The mapping between |l9| ) and the loop model 
follows from a straightforward graphical expansion which 

is described in 133- 


Let us consider the continuum description of (19). To 


begin with, take a configuration in which z is slowly vary¬ 
ing. Each term in the product over nodes is then close 
to one, and we may obtain a continuum sigma model 
Lagrangian by a derivative expansion. In 3D, the only 
term with two derivatives allowed by global, gauge and 
lattice symmetries is the standard sigma model kinetic 
term. Let us focus on the n = 2 case, and parameterise 
CP 1 (which is simply the sphere) using the Neel vector 


M a = Z^CTaZ, 


a = x, y, z, 


( 20 ) 


instead of the gauge-redundant field 2 . Then 

(iV 2 = 1). (21) 


C a = y(ViV) 2 , 


A crude way to estimate a bare value of K is to calculate 
the Boltzmann weight in (19) for a spin configuration 


with a uniform twist, giving K = 1/16. (For general n, 
L a may be written in terms of the matrix Q = zz' — 1/n.) 

The Lagrangian obtained by the derivative expansion 
can fail to capture the true scaling behaviour in two ways. 
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It fails in a trivial way when N varies strongly at a node. 
This will of course be the case in the lattice model, and 
leads to an order one renormalisation of the stiffness. 

More importantly, the phase of z can vary rapidly even 
if N is slowly varying. Nodes where N is approximately 
constant but where this phase varies abruptly contribute 
imaginary terms to the action. For smooth configura¬ 
tions of trivial topology, these phases cancel. However, in 
the presence of hedgehog defects, there remain nontrivial 
phases that are missed by the derivative expansion. 

This is because in a configuration with a hedgehog it 
is impossible to find a gauge in which z is everywhere 
slowly varying, even far from the hedgehog core. This 
follows from the fact that topological flux density 
which when integrated over a closed surface gives the 
signed number of hedgehogs inside, is a total derivative 
when written in terms of 2 ;: B M ^{z^\/^z). 

If z were continuous, integrating the topological density 
over a large sphere would give zero. Therefore z must 
be discontinuous somewhere on the sphere if the sphere 
encloses a hedgehog. 

A simple calculation is required to determine what ef¬ 
fect the imaginary terms have on the weight for a con¬ 
figuration with a hedgehog. We do this calculation in 
Appendix [B] For specificity, we take the hedgehog to be 
centred on a site of C\ or a site of C 2 — for example at the 
centre of the cube in Fig. [T] These locations form a bcc 
lattice, with four sublattices. We find that the weight 
of a configuration with a hedgehog acquires a fugacity 
proportional to 

1 , i, —1 or — i, ( 22 ) 

depending on which sublattice it sits on. (More pre¬ 
cisely, only the relative phase between different locations 
is meaningful [51].) 

This also generalises immediately to larger n. The re¬ 
sult matches nicely what is found for the 2D quantum 
Heisenberg model [55H551 . where the fugacity for instan- 
tons hedgehogs in spacetime — takes the same set of 
values as above depending on which of four sublattices of 
the square lattice the instanton occurs on. 

By symmetry, we infer that the coarse-grained hedge¬ 
hog fugacity vanishes. That is, it vanishes as a result of 
phase cancellation between configurations in which the 
hedgehog is centred on nearby sites on different sublat¬ 
tices. Thus the arguments of Refs. CHI apply, giving the 
NCCP 1 model (Eq. [lj as the continuum description. 

An alternative argument for the NCCP 1 description 
of square lattice spin- 1/2 antiferromagnets was given in 
Ref. [55], focusing on the VBS order parameter (p and its 
vortex defects rather than the Neel order parameter N 
and its hedgehog defects. The key point of this alterna¬ 
tive argument is that a vortex in the VBS order parame¬ 
ter carries a single unpaired spin at its centre. In space- 
time, this corresponds to an extended spinon worldline 
running along the vortex core. This has a direct interpre¬ 
tation in the loop model: a vortex line in the node order 
parameter (p has a single extended loop running along its 


core. We have confirmed this explicitly by constructing 
such configurations. 

In previous work we have considered transitions in a 
different version of the loop model which does not pre¬ 
serve the full lattice symmetry (36) [37]. The transitions in 
that less-symmetric model are described by the compact 
CP " -1 model, unlike the present loop model whose tran¬ 
sition is described by NCCP" -1 . This is because break¬ 
ing lattice symmetry spoils the cancellation between the 
values in Eq. [22] leaving a nonzero hedgehog fugacity. 
This allows the standard critical behaviour of the com¬ 
pact CP " -1 model (i.e. of the usual 0(3) model when 
n = 2). We note that in the case n = 3, the compact 
CP 2 model appears to show an interesting continuous 
transition which is naively forbidden by Landau theory 
[36l 37] [56] ■ An RG explanation for why the expectation 
from Landau theory breaks down in this case was given 
in Ref. [57] . 


VI. RG FLOWS IN THE NCCP" -1 MODEL 

(n-COMPONENT ABELIAN HIGGS MODEL) 

In this section we make a conjecture for the topology of 
the RG flows in the NCCP " -1 model, the ?r-component 
generalisation of Eq. [I] with z = (zi,..., z n ): 

C = |(V — iA)z\ 2 + k(V x A) 2 + /i|z | 2 + A|z| 4 . (23) 

This is also known as the n-component Abelian Higgs 
model. We will treat both n and the spatial dimension d 
(between 2 and 4) as continuously varying. 

Scaling violations are seen in a wide variety of differ¬ 
ent lattice models that are related to this field theory (at 
n = 2 ) and persist to very large lengthscales [57] [53], so 
we believe a plausible explanation for them should appeal 
to universal physics of the NCCP " -1 model and not to 
accidental features of specific Hamiltonians. Results for 
SU(3) and SU(4)-symmetric models (n = 3, 4 [10, 27]) 
suggest that a satisfactory explanation should also ac¬ 
count for scaling variations across a range of n. 

Fig. [L9] shows the basic topology of the RG flows we 
find. This is a sheet of RG fixed points projected on the 
space of n, d, and a scaling variable A. (Close to 4D, A 
is the quartic coupling in Eq. |T] but in lower dimensions 
one cannot make this identification.) The RG flow is 
parallel/antiparallel to the A axis, since n and d do not 
flow. A is irrelevant on the critical sheet and relevant 
on the tricritical sheet. The strongly relevant coupling 
which drives the transition (i.e. the mass) is not shown, 
since we consider the theory at the critical value. For 
any fixed value of n between zero and n *4 ~ 183, the 
critical point exists so long as we are sufficiently close to 
two dimensions. When d is increased, the critical point 
disappears at a universal value d*(n), by merging with 
the tricritical point. 

To begin with, consider three limits in which the 
NCCP " -1 model is solvable. First, it is tractable by 
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FIG. 19. Figure showing topology of the sheet of RG fixed 
points for dimensionalities around three (the limiting cases 
d = 2 and d = 4 are not shown). 


saddle-point at large n, where it yields a nontrivial crit¬ 
ical point for 2 < d < 4. This critical point describes a 
direct transition between a Higgs phase, where 2 is con¬ 
densed and SU(n) symmetry is broken, and a Coulomb 
phase where SU(?z) symmetry is unbroken and the gauge 
field A is massless. 

The field theory is also tractable in a 4 — e expan¬ 
sion [591 . For infinitesimal e, a weak-coupling critical 
point exists only if n is greater than or equal to a value 
which we denote n* 4. This value is quite large, 77* 4 ~ 183. 
In fact in the regime n > 77*4 where the critical point ex¬ 
ists, the 4 — e expansion also yields a tricritical point 
at a smaller value of the quartic coupling A. As n ap¬ 
proaches n *4 from above, these two fixed points approach 
each other, and they annihilate when 77 reaches ?7*4. For 
n < 77*4 there is no nontrivial fixed point: the theory is 
expected to flow to a discontinuity fixed point at large 
negative A representing a first order transition. 

Finally, the NCCP n_1 model can be studied in 2 + 
e dimensions by switching from a soft-spin formulation 
to a nonlinear sigma model [60l 161] . In this regime a 
continuous phase transition is found for all values of n 
greater than zero. (The ‘replica-like’ regime n < 1 is 
meaningful and describes certain classical loop models 
[36l 62, 63].) In the 2 + e approach it does not matter 
whether the nonlinear sigma model is formulated with a 
dynamical U(l) gauge field, or as a pure nonlinear sigma 
model with target space CP n_1 : the two formulations 
give identical results (60] . (When the dynamical gauge 
field is included its coupling flows to infinity, so it can be 
integrated out leaving the usual CP ra_1 nonlinear sigma 
model.) 

A crucial point, made in Ref. [3T], is that the fixed 
point found in the large 77 approach is the same as that 
found in both the 2 + e and 4 — e expansions. This can be 


seen by comparing the results for the critical exponents in 
the regions of overlap of the expansions. Viewing n and 
d as continuous variables, there is therefore a continuous 
family of fixed points in a region of the (n, d) plane for 2 < 
d < 4 and sufficiently large 77 ED- This region is defined 
by n > n*(d), where n»(d) is the d-dependent value of 
n at which the fixed point disappears. From the 4 — e 
approach we know that the limiting behaviour as d —> 4 is 
n*(d) —» n *4 ~ 183, and from the 2+e expansion we know 
that 77* (d) —)• 0 as d —>■ 2. Using the sigma model we may 
also argue that the slope of 71* (d) is finite as d —>■ 2 (see 
endnote [M]). The value of 77* (3) is no t know n, but it 
is possible to show that 71 *(3) > 1 (Sec. VI 3). We will 
assume that 77* (d) increases monotonically with d. 

We can go further by noting that in the 4 —e expansion 
the way in which the nontrivial fixed point disappears at 
77* (d) is by annihilation with a tricritical point. On the 
basis of continuity we expect that the mechanism for the 
disappearance of the fixed point at 77 * (d) is the same for 
all d. This leads directly to the topology in Fig. HI It 
will be convenient to denote the value of d where the 
merging of the critical and tricritical points happens, for 
a given 71, by d*(r 7 ). 


1. RG flows close to merging line 


Since 77*4 ~ 183 is relatively large, the average slope of 
the line d = d^(n) is small. It is therefore possible that 
there is a broad range of 71 -values where the line lies close 
to three dimensions, i.e. where |d*(n) — 3| is small. So it 
is worth studying the RG flows in this regime. 

On the line, A is marginal. After rescaling and shifting 
A by a constant, its RG equation is: 


dA 

d In L 


-A 2 . 


(24) 


Moving slightly away from the line, the RG equation be¬ 
comes, to lowest order in d*( 77 ) — d, 


-j-^r ~ 0(77) (d* (77) - d) - A 2 ( 25 ) 


with an unknown but universal positive constant a(n). 
This equation encapsulates the fact that when d > d* (77) 
there is no fixed point, and when d < d* (77) both a critical 
and a tricritical point exist, at A = ± \Ja(n ) (d* in) — d). 

Fixing now on d = 3 , we define the universal quantity 
A(77) = a(n)(d*(n)— 3 ), which is zero at 77 = 77 *( 3 ) (where 
the critical point disappears in 3D) and small over the 
range of 77 where the merging line lies close to d = 3 : 


dA 

dlnL 


A (77) — A 2 . 


(26) 


When 77 < 77* (3), A(n) is negative and there is no fixed 
point in 3D. Instead the RG flows go off to large negative 
A, suggesting a first order transition. However if we are 
close to the line, so |A(n)| is small, the RG flow becomes 
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very slow. Integrating Eq. [26] shows that this implies an 
exponentially large correlation length at this transition: 

£ ~ exp (tt/vIA( n)|) (27) 


3. Bound on n*(3) 

The value of n*(3) is not known, but we can argue that 
n*(3) > 1. (29) 


with a nonuniversal prefactor. (A similar phenomenon 
occurs in the 2D Q-state Potts model for Q > 4 [551 .) 

For n > n*(3) there is a conventional critical point at 
A = yf A(n). But if we are close to the line the lead¬ 
ing irrelevant RG eigenvalue, y- lrr , at this critical point is 
small, implying large corrections to scaling: from Eq. 25 


y irr ~ -2a/A( n). 


(28) 


2. Interpretation 


At first sight the topology we have found for the RG 
flows suggests two possible explanations for the scaling 
violations. 

First, if we (speculatively) assume that ro*(3) > 4, 
but that the merging line lies close to d = 3 over the 
range 2 < n < 4, then we obtain an anomalously weak 
first order transition for n = 2,3,4. In this scenario 
there is pseudo-critical behaviour (with drifting expo¬ 
nents EZI) up to the exponentially large lengthscale of 
Eq. [27j thanks to the ‘nearby’ fixed point at slightly 
smaller spatial dimension. The virtues of this scenario 
are that it appeals to universal features of the RG flow, 
so may explain why numerous different lattice models see 
very similar scaling violations, and that it can produce 
scaling violations over a range of n. 

(Note: if |A| is very small, there exists a range of sizes 
where A appears to be a conventional marginally irrele¬ 
vant variable. However this range ends at a size L* which 
is parametrically smaller than £ [681 when £ is large. The 
basic point is that the stages of the RG flow with A > 0 
and with A < 0 take roughly equal amounts of RG time, 
but the latter corresponds to a vastly larger lengthscale 
because of the logarithmic relationship between length 
and RG time.) 

Second, we might try to explain the scaling violations 
differently by postulating that in 3D the values n = 2, 3,4 
lie just below the merging line (n*(3) < 2) so that A(ra) 
is small and positive for n = 2,3,4. This would give a 
true critical point with large (but conventional) scaling 
corrections due to a small irrelevant exponent y-„ r . How¬ 
ever, our numerical results strongly indicate that a small 
y 1TI is not sufficient, on its own, to explain what we see. 

We emphasise that since our argument fixes the topol¬ 
ogy of the RG flows but not the numerical value of n* ( d ), 
the scenario above for a weak first order transition is 
speculative. In Sec. |VII| we will discuss other possible 
conjectures. 


This may look surprising at first glance, since the single¬ 
component Abelian Higgs model, n = 1, certainly has 
a continuous phase transition in 3D. This transition is 


means that this ‘inverted XY’ phase transition does not 
lie on the sheet we are considering: it is not analytically 
connected to the deconfined critical point at large 11 m- 

Formally one can see this as follows. If the inverted 
XY transition did lie on the critical sheet of Fig. [19J 
we could describe it by setting n = 1 in the 2 + e ex¬ 
pansion of the CP™ -1 nonlinear sigma model. But this 
is evidently not the case. The inverted XY transition 
is a conventional thermodynamic phase transition with 
nontrivial signatures in the free energy. In contrast, the 
sigma model at n = 1 is a replica-like theory, in which 
the number of degrees of freedom becomes zero and the 
free energy vanishes identically. The same reasoning im¬ 
plies that n*(3) is strictly greater than one. Otherwise 
the n = 1 model would have a Higgs transition with an 
unphysical replica-like continuum description. 

Instead of being connected to the critical points in 
Fig. [l9j the inverted XY transition is the n = 1 limit 
of a much simpler transition - namely that which (for 
n > 1) separates the Higgs phase, where z is condensed, 
from a pair condensed phase where the bilinear zz' is 
condensed but z itself is not. (This phase appears for 
appropriate couplings [15][71].) SU(n) symmetry is bro¬ 
ken on both sides of this transition, so it is not like the 
critical points in Fig. [19] One can check that the tran¬ 
sition is in the inverted XY universality class for all n 
(as the interactions between the critical sector and the 
Goldstone modes of the broken SU(n) symmetry are ir¬ 
relevant), so the n dependence of the critical behaviour 

is trivial IE- 


related by duality to that of the XY model 


Eq. 29 


4- ‘Failure’ of the 2 + e expansion of the 0(3) model 

The CP 1 nonlinear sigma model is the 0(3) nonlinear 
sigma model by another name, as the target space is the 
sphere (CP 1 = S 2 ). Therefore the topology of the flows 
in Fig. [19] confirms that the standard 2 + e expansion 
of the 0(3) sigma model does not describe the Wilson 
Fisher critical point of the 3D 0(3) model. Instead, set¬ 
ting e = 1 in this expansion describes the SU(2) decon¬ 
fined critical point, if it exists (or nothing at all if the 
critical point vanishes at a d*(2) below 3). The applica¬ 
bility of the 2 + e expansion to the abelian Higgs model 
was also argued for in Ref. m- 

Although this is contrary to what is often assumed, it 
should not be surprising in the light of knowledge about 
hedeghogs in the 3D 0(3) model g] 173] , Suppressing 
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these topological defects has been convincingly argued 
to change the critical behaviour [4] . It is natural that 
the 2 + e expansion, which considers only spin waves, 
describes the behaviour in the absence of hedgehogs m- 

The conclusion is also supported by an RG approach 
to the O(M) model that employs a double expansion in 
(d — 2) and (M — 2) [32] . This shows that the standard 
2 + e expansion fails to capture the critical behaviour of 
the O(M) model when M is smaller than a d-dependent 
critical value M c . To first order in d— 2, the relationship is 
(M c — 2) ~ — 2), which gives M c ~ 4.5 in 3D. While 

higher order corrections in (d — 2) could be significant, 
this result suggests M c > 3 [TH? . 

Ref. =22] also notes the poor agreement between the 
2 + e exponents in 3D and the exponents of the 0(3) 
model in 3D. For example, v is equal to v = 1/2 at order 


2 , or to 2/5 at order 


EEj. In the light of the 


preceding, we should instead apply these exponents to 
the SU(2) deconfined transition. Indeed the e-expansion 
values v = 1/2 and v = 2/5 are remarkably close to our 
best estimates of v (see Sec. IVE). 


VII. INTERPRETATION AND CONCLUSIONS 


We have argued that the scaling violations at the de¬ 
confined critical point are too severe to be explained as 
corrections to scaling from a weakly or marginally irrel¬ 
evant scaling variable, and have sharpened the possible 
alternatives. 

Some of our numerical results suggest that estimates 
of the exponents v, ?7Neei and t?vbs may be better defined 
if we avoid using finite-size scaling to obtain them (see in 
particular Secs. IV A 1 IV E) — although of course aban¬ 


doning finite size scaling restricts us to lengthscales much 
smaller than the system size. The most radical conjecture 
would therefore be to attribute the scaling violations to a 
dangerously irrelevant variable (DIV) which leaves criti¬ 
cal behaviour intact but modifies finite-size scaling (see 
also Ref. m • In the simplest picture, the role of this 
DIV would be to cut off fluctuations of some zero-mode(s) 
of the fields that are unbounded in the pure fixed-point 
theory. The main examples we know of this phenomenon 
are in free theories (such as /> 4 theory above 4D j20j or the 
quantum Lifshitz theory CZD; another type of example 
is in theories that are dual to free theories [78]. Further 
work is required to determine whether it is a plausible 
possibility in an interacting theory such as Eq. [T] 

The alternative possibility of an anomalously weak first 
order transition has been discussed in detail in Sec. eh 
W e have seen that in principle there is a mechanism by 
which a very large correlation length can appear without 
the need for fine-tuning of the Hamiltonian, and that 
in this scenario there would be pseudo-critical behaviour 
over a large range of scales with (for example) drifting 
critical exponents. Such a possibility is hard to exclude 
numerically. We note however that we do not see the 
usual signs of an incipient first order transition, despite 


studying much larger scales than the early simulations 
used to argue for first-order behaviour in the J-Q model 
M- The first-order scenario would also leave the good 
scaling of, for example, the derivatives of the correlators 
(Sec. IV A 1) a mystery. 

Intriguing questions therefore remain for the future. 
The loop model is an ideal platform for further work on 
the deconfined transition, since it provides an intuitive 
geometrical picture and since isotropy in three dimen¬ 
sions is a convenient feature. It would also be interesting 
to perform simulations at other values of n (which in the 
formulation just after Eq. [3] need not be an integer) in 


order to probe the scenario of Sec. VI 


Various modifications to the loop model are possible. 
For example one may allow a third node configuration 
(see Fig. [2]) in which the two incoming links are joined 
and the two outgoing links joined, so that the loops are 
no longer consistently oriented. The symmetry is then 
broken from SU(n) to SO(n), and the relevant field the¬ 
ories are RP” -1 models [Ml E3 m 023 which can show a 


Z 2 spin liquid phase [80] [81]. One may also study the ef¬ 
fects of various anisotropies and symmetry-breaking per¬ 
turbations. Finally, the loop model generalises to four 
dimensions, where it may be a useful tool to search for 
new types of critical behaviour arising in 3+ID quantum 
magnets. 
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Appendix A: Methods 

The numerical procedure is as follows. An initial state 
is constructed by choosing at random one of the two pos¬ 
sible configurations of a node with equal probabilities. 
Loops are then formed by following turning instructions 
at each node. We generate the fugacity n from a sum on 
loop colours, assuming n is integer. So a colour is asso¬ 
ciated with each loop chosen with equal probability from 
n alternatives. 

Subsequent states are generated using three kinds 
of parallelized Monte Carlo moves to ensure that at 
equilibrium configurations are distributed according to 
the partition function Eq. [4] The first one up¬ 
dates the state of the nodes using a checkerboard¬ 
like algorithm. Each node, if its four links have 
the same colour, changes its state with probability 
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minjexp — 2Ja(i) Ylj nn a ti) , lj- The updates of all 
the nodes in the lattice are done in three stages. In the L 
lattice there are two sublattices A and B. Each sublattice 
is tripartite and the three sublattices are simple cubic. At 
each time one sublattice of A and one sublattice of B are 
updated. The second type of Monte Carlo move chooses 
a link at random and changes the colour of all the links 
of the associated loop to a different colour, chosen with 
uniform probability from the n — 1 possibilities. This 
move is also parallelized by letting each thread choose a 
link and change the color of the loop if it is not already 
visited. The third type of move is to re-colour all loops in 
the system, with the new colours selected independently 
and at random for each loop. It is designed to ensure 
that the colours of short loops equilibrate effciently. 


A number of these moves are combined to form a com¬ 
posite update which we term a Monte Carlo sweep. The 
ingredients in a single sweep are as follows. First we it¬ 
erate 20 times a sequence in which moves of the first two 
types are intercalated and repeated for each of the three 
sublattices. Then we apply the third type of move. Mea¬ 
surements are performed once every Monte Carlo sweep. 
The autocorrelation function of the energy is used to es¬ 
timate a correlation time. The blocking and bootstrap 
methods [82] are used to estimate errors. 

As an interpolation scheme we have used Ferrenberg’s 
multiple histogram method [53] to obtain a continuous 
set of values as a function of J of the different quantities. 
This method is also employed whenever a derivative has 
to be calculated and a maximum or minimum has to be 
obtained. Errors related to this technique are calculated 
using the bootstrap method. 

We consider system sizes of up to 3.9 x 10 s links or 
lateral size L = 640, with extensive results for L < 512. 
The minimum number of Monte Carlo sweeps used is 10 5 
for any J and L , and increases with decreasing L. 


Appendix B: Calculation of hedgehog fugacity 


->n—1 


Regarded as a lattice magnet for classical CP r " ' spins, 
the partition function (191 has the peculiar feature that 
although it is both local and gauge invariant, it is not 
simply expressed in terms of the gauge invariant quan¬ 


tity N (for n = 2) or Q (for larger n). The consequence 
of this, as noted in Sec. [Vj is that the sigma model ac¬ 
tion arising from a derivative expansion may need to be 
supplemented by purely imaginary terms from nodes at 
which the phase of z changes abruptly. In the presence 
of hedgehogs such nodes are inevitable, even far from 
the hedgehog core. As a result the hedgehog fugacity 
acquires a spatially varying phase. 

For simplicity we take the hedgehogs to sit at the centre 
of a cube of C\ or Ci (Eq. [2]), for example at the centre 
of the cube in Fig. [l] (left). These locations form a bcc 
lattice. We take the origin at one such bcc site and the 
coordinate axes parallel to the links. 


Focussing on the case n = 2 (the generalisation to 
larger n is immediate), let us first consider the representa¬ 
tive configuration in which the hedgehog is centred at the 
origin and the Neel vector N (defined in Sec. [v]) points 
in the radial direction. In polar coordinates, 0 < 9 < n, 
0 < ip < 2-7T, this is 

N = (sin 9 cosy), sin 9 sin <p, cos 0), (Bl) 


where the coordinates of a link are those of its midpoint. 
To write this in terms of z we must pick a gauge. It is 
convenient to choose one in which the Boltzmann weight 


-s n 


1 - 
2 . 


0 l z i)i z l'Zi') + {zlzi'){zl,Zi)] (B2) 


is approximately equal to one for as many nodes as possi¬ 
ble. For the links with positive and negative 2 coordinate 
(9 < 7t/ 2 and 9 > 7r/2) we take, respectively, 

z = (cos 9/2, e llp sin 9/2), z = {e~ q ' v cos 9/2, sin 9/2). 

(B3) 

We see from Fig. [l] that there are also links in the equa¬ 
torial plane. For these we take 


z = -^{e~ ilfi/2 ,e iip/2 ). 

v 2 


(B4) 


Now consider e -15 " 0110 . In fact it will suffice to consider 
only nodes far from the core, and to treat z as constant 
except for the discontinuities in our gauge choice: other 
contributions to the action are either included in the spa¬ 
tially independent amplitude of the hedgehog fugacity, or 
are already captured by the naive derivative expansion. 

With the above gauge choice, the nodes at which z 
varies abruptly all lie in the equatorial plane, and have 
two of their links within this plane, one above it, and 
one below. Fi g. [I] shows four such nodes (all black). 
From Eqs.|B2|- B4 we find that most of these nodes have 


e -S n ode ~ x, despite the variation in the phase of z. How¬ 
ever there is a string of nodes along the positive x axis 
(ip = 0, 9 = 7t/ 2) each of which contributes a minus sign 
to the Boltzmann weight (i.e. e~ SllodB ~ —1 for these 
nodes). If we translate the core of the hedgehog by the 
vector (2,0,0), while keeping the configuration far from 
the core fixed, we change number of nodes on this string 
by one. Therefore this translation changes the sign of the 
Boltzmann weight. 

Let the phase term in the hedgehog fugacity be denoted 


exp (i rj{r)), 


(B5) 


where the spatial vector r lies on a bcc site. (An ‘anti¬ 
hedgehog’ of negative topo logic al charge has phase factor 


e as we see from Eq. B2 and the fact that complex 


conjugating z exchanges hedgehogs and antihedgehogs.) 
The phase ?y(r) is defined only up to a constant: for exam¬ 
ple in a closed system there are equal numbers of hedge¬ 
hogs and antihedgehogs, so the constant part of r] drops 
out. 
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It may be seen from Fig. [T] that the translational sym¬ 
metry between bcc sites is not spoiled by the link ori¬ 
entations. Using this, we may argue that rj(r) is of the 
form rj(r) = k.r for some momentum k. By the above 
calculation, e lri U) = — e z? h r + 2x ), where x = (1,0,0). By 
symmetry, we have similar results for translations in the 


y and z directions. This is enough to fix k up to a sign: 

* = ±£( 1,1,1)- (B6) 

One of these signs applies to the hedgehog and one to the 
anti-hedgehog. We have not fixed which is which, but it 
does not matter. 

With this k, the hedgehog fugacity takes four distinct 
values on the four sublattices of the bcc lattice, propor¬ 
tional to ±1 and ±i. This is the result quoted in Sec. [V| 
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